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Abstract 

In three-dimensional O(A^) models, we investigate the low-momentum be- 
havior of the two-point Green's function G{x) in the critical region of the 
symmetric phase. We consider physical systems whose criticality is charac- 
terized by a a rotationally-invariant fixed point. 

In non rotationally-invariant physical systems with 0(A^)-invariant inter- 
actions, the vanishing of anisotropy in approaching the rotationally-invariant 
fixed point is described by a critical exponent which is universal and is 
related to the leading irrelevant operator breaking rotational invariance. At 
N = oo one finds p = 2. 1/N expansion and strong-coupling calculations 
show that, for all values of > 0, /? ~ 2. 

Non-Gaussian corrections to the universal low-momentum behavior of 
G{x) are evaluated, and found to be very small. 

PACS numbers: 05.70.Jk, 64.60.Fr, 75.10.Hk, 75.40.Cx 
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Three-dimensional 0{N) models describe many important critical phenomena in nature. 
We just mention that the case N = 3 describes the critical properties of ferromagnetic 
materials. The case = 2 is related to the helium superfiuid transition. The case = 
1 (i.e., Ising-like systems) describes liquid- vapour transitions in classical fluids or critical 
binary fluids. Finally, the limit — > is related to dilute polymers. 

The critical behavior of the two-point correlation function G{x) of the order parameter is 
related to the phenomenon of critical scattering observed in many experiments, e.g. neutron 
scattering in ferromagnetic materials, light and X-rays in liquid-gas systems. In Born's 
approximation, the cross section Tfi for incoming particles (i.e. neutrons or photons) of 
momentum pi and final outgoing momentum pj is proportional to the component k = Pf—Pi 
of the Fourier transform of G{x) 

Tfi(xG{k), k = pf-pi. (1) 

As a consequence of the critical behavior of the two-point function G{x) at Tc, which is 

Gik) - (2) 

the cross section for A; ^ (forward scattering) is observed to diverge as T — T^. When 
strictly at criticality, the relation (|^) holds at all momentum scales -C A, where A is a 
generic cut-off related to the microscopic structure of the statystical system. In the vicinity 
of the critical point, where the relevant correlation length ^ is large but finite, the behavior 
(I) occurs for A > A; > 1/^. 

We will specifically consider systems with an 0(A^)-invariant Hamiltonian in the symmet- 
ric phase, i.e. where the O(A^) symmetry is unbroken. Furthermore, we will only consider 
systems with a rotationally-symmetric fixed point. Interesting members of this class are 
systems defined on highly symmetric lattices, i.e. Bravais or two-point base lattices with a 
tetrahedral or larger discrete rotational symmetry. 

Let us consider the second moment correlation length 

^•2 _ 1 X^x |xp(j(x) , , 

and the mass scale Mq = ^/^g- In the critical region of the symmetric phase and at low 
momentum, experiments show that G{x) is well approximated by a Gaussian behavior, i.e. 

^ ~ 1 H (41 

G{k) ~ Ml ^> 

Our aim is to estimate the deviations from Eq. @) in the critical region of the symmetric 
phase, i.e. for < T/T^ — 1 <C 1, and in the low-momentum regime k"^ ^ Mq. We focus on 
two quite different sources of deviations: 

(i) Scaling corrections to Eq. (^), depending on the ratio A;^/M^, and reflecting the 
no n- Gaussian nature of the fixed point. 

(ii) Non rotationally-invariant scaling violations, reflecting a microscopic anisotropy in 
the space distribution of the spins (assuming that no anisotropy is generated by their in- 
teraction). This phenomenon may be relevant, for example, in the study of ferromagnetic 
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materials, where the atoms he on the sites of a lattice, and anisotropy may be observed in 
neutron scattering experiments. In these systems anisotropy vanishes in the critical limit, 
and G{x) approaches a rotationally-invariant form. 

We mention that the effects of a breakdown of O(A^) symmetry in the interactions have 
been widely considered in the literature 

Several approaches have been considered in order to study the critical behavior of the 
two-point function 

G{x) = {s{x)-m)- (5) 

In lattice O(A^) non-linear a models, we have calculated the strong-coupling expansion of 
G{x) up to 15th order on the cubic lattice and 21st order on the diamond lattice within 
the corresponding nearest-neighbor formulations, 

Sl = -N/3Y.s^^-s,^, (6) 

links 

where xi,Xr indicate the sites at the ends of each link. Our conclusions are integrated and 
supported by results obtained from a. 1/N expansion of lattice O(A^) models. Furthermore, 
we have analyzed the first nontrivial term of the e-expansion and of the gf-expansion (i.e. 
expansion in the coupling at fixed dimensions = 3) of the two-point function within the 
corresponding 0^ formulation of O(A^) models. 

Our results substantially confirm experimental observations that scaling corrections to 
the Gaussian behavior (^) are very small. This point has been already discussed in the liter- 
ature, essentially by e-expansion and 1/A^ expansion |Q in the continuum formulation of 
0{N) models. We present here rather accurate determinations of such corrections obtained 
by a strong-coupling analysis of G{x). 

We will show that the anisotropy of G{x), for the class of systems we are considering, 
vanishes at the rotationally-invariant fixed point, with a behavior governed by a universal 
critical exponent p: non-spherical moments (i.e. those which vanish when calculated on 
spherical functions) of G{x) are depressed with respect to spherical moments carrying the 
same naive physical dimensions by a factor From a field-theoretical point of view, 

anisotropy in space is due to non rotationally-invariant (but O(A^) symmetric) irrelevant 
operators in the effective Hamiltonian, whose presence depends essentially on the symmetries 
of the physical system, or of the lattice formulation. Another possible origin is the coupling 
of the system to an external (anisotropic) source. The exponent p is related to the critical 
effective dimension of the leading irrelevant operator breaking rotational invariance. On 
cubic-like lattices the leading operator has canonical dimension d + 2. In the large- limit, 
where the canonical dimensions determine the scaling properties, one then finds p = 2. 1/A^ 
expansion and strong-coupling calculations show that, for all values of > 0, p remains 
close to its canonical value. One can also show that for the two-dimensional Ising model 
the exact result p = 2 holds. Moreover, in the 1/A^ expansion of two-dimensional models, 
one finds that, to 0(1/A^), the multiplicative logarithmic correction to the behavior ~ 
is absent. 

The technical details of our study will be reported in a separate extended paper. Here 
we only describe the general features of our analysis and the main numerical results. 
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For definiteness let us consider the cubic lattice version of the models. We parametrize 
the two-point spin-spin function by a multipole expansion in the form 



p-^G-\k,MG) = Y.92i{y,MG)Q2i{k), (7) 

where y = k'^/MQ, and Q2i{k) are homogeneous functions of momenta of degree 21 which 
are invariant under the symmetries of the lattice. Their expressions can be obtained from 
the fully symmetric traceless tensors of rank 21, T2i^"''^^'{k), by considering all the cubic- 
invariant combinations [^. Odd rank terms are absent in the expansion (|^) because of 
the parity symmetry. Moreover, there is no rank-two term, due to the discrete rotational 
symmetry of the lattice. The first nontrivial function is 

Q4(fc)=E<-|(E^;)'- (8) 

Q4:{k) corresponds to the leading irrelevant operator that breaks rotational invariance in the 
low-momentum expansion of the Hamiltonian. In the continuum notation this operator has 
the form 

0^{x) = s{x) ■ Q4{d)s{x), (9) 

which has canonical dimension five in d = 3. 

The scaling limit of Eq. (|^) corresponds to taking Mq while keeping the ratio k/Mc 
finite. Hence we can parametrize the low- momentum behavior of G~^{k, Mq) in the critical 
region by 

p-'G-\k,MG) = Z'%{y)M^G + - 

+ ZMk/MG)g,{y)M^ + 0{Ml;), (10) 

where we have dropped rotationally-invariant 0{Mq) terms which we are not interested in, 
and we have introduced the quantities 

= — ' ^4 = c/4(0,Mg), (11) 

G 

which, in the limit Mq 0, absorb all non-analytical dependence on Mq of the correspond- 
ing terms. The functions go{y) and g^i^y) are universal, i.e., they do not depend on the 
specific form of the lattice hamiltonian. They possess a regular expansion around y = 0: 



goiy) = l + y + J2c^y\ (12) 

OD 

g,{y) = l + J2d,y\ (13) 



i=l 



In the limit N ^ oo, the models are strictly Gaussian and therefore all the coefficients 
c,- and di are zero, 
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We consider the spherical moments 



m2j = J2\^\^^G{x), (14) 

X 

and the leading non-spherical moments 

94,i = E \x\^'Q4{x)G{x) (15) 

X 

which vanish if G{x) is rotationally invariant. The critical exponent p, describing the van- 
ishing of anisotropy, and the coefficients q and di of the low-momentum expansion of go{y) 
and g^iy), can be determined by studying appropriate combinations of the above moments 
in the critical limit. In the critical region 

~ ^, (16) 

where g2,m and m4+2m have the same naive physical dimensions. 

In the case of non-Gaussian fixed points, like those corresponding to the theory at finite 
N, O^lx) develops an anomalous dimension cr, which causes a departure from the Gaussian 
value of p, p = 2 + a. a can be extracted by evaluating the ratio 

Mi g„(0,Mc) ° ^ ' 

In turn this combination is easily estimated by taking the moment ratio g4,o/'^2 ~ Mq. We 
have calculated the first non-trivial term of the 1/N, e and (^-expansions of a: 

- = JL + 0(±), (18) 



21tt'^N 



a = g*^ T 7^ + o*^ (20) 

^ 25515 {N + Sy ^ ^ ^ 

{g* is the fixed-point value of the rescaled coupling: g = g{N + 8)/4:8'k [^). Strong-coupling 
estimates of a have been obtained by analyzing and comparing the available strong-coupling 
series of q^^Q and m2 on both cubic and diamond lattices. Universality between cubic and 
diamond lattice is substantially verified, although the analysis on the diamond lattice leads 
to less precise results. In Table | we report the strong-coupling estimates of a, and for 
comparison its estimates from Eqs. (0), (|T^ and (PDj). For O(A^) models, the values of a 
are very small in the whole range > 0, thus indicating an essentially Gaussian behavior 
of this critical exponent. The comparison of the estimates of a from the various approaches 
we have considered is satisfactory. 

The low-momentum critical behavior of the two-point Green's function, i.e. 
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(21) 



has been already investigated within an e-expansion of the 0^ formulations of O(A^) models 
(up to 0(e^)) 0, and in the expansion (up to 0(l/iV)) The computations of the 
0(l/iV^) contribution to go{y) is in progress. Moreover rough estimates of such corrections 
have been presented for = 1 and iV = 3 by analyzing strong-coupling expansion (up to 
10th order) within various lattice formulations of the theory . We have reconsidered the 
problem of determining the non-Gaussian corrections to go{y) in the low- momentum regime, 
and those of giiy), especially by a strong-coupling analysis. Furthermore we have calculated 
the first nontrivial contributions to ^4(1/) in the 1/N, e and (7-expansions, and of go{y) in the 
(7-expansion, and the corresponding coefficients of the expansion around y = 0. 

Our strong-coupling analysis leads to a satisfactory precision which, in the case of the 
coefficient q, considerably improves earlier calculations. This has been achieved essentially 
for two reasons: longer strong-coupling series are available, and, even more important, 
improved estimators have been employed. We indeed took special care in the choice of 
estimators for the physical quantities c, and dt. This is very important from a practical 
point view: better estimators can greatly improve the stability of the extrapolation to the 
critical point. Our search for optimal estimators was guided by the knowledge of the large- 
limit. We chose estimators which are perfect for = 00, i.e. do not present off-critical 
corrections to their critical value q = cij = in the symmetric phase. In particular in the 
case of C2 and di and on the cubic lattice, we used the following estimators: 

M^ + Im^, (22) 




120mo ^ 20 

^4,1 .^2 , 1 



M^ + —M^] (23) 

44^4,0 ^ U ^) ^ ^ 

respectively. In Table |T| we report our strong-coupling estimates for C2, C3 and di on both 
cubic and diamond lattice, obtained by evaluating at Pc appropriate approximants (such as 
Fade and first order integral approximants) of the strong-coupling series of the corresponding 
estimators. One may notice that universality between cubic and diamond lattice is always 
confirmed. At large- A^ the strong-coupling estimates of C2, C3 and di compare very well 
with the corresponding 0(1/A^) formulas: C2 = -0.004449/A^ + C3 = 0.000134/A^ -t- 
etc... 0, di = -0.0020647/A^ + ^2 = 0.0000738/A^ + etc.... The comparison with the 
estimates from other approaches, which are reported in Table 0, is good. 
All calculations agree in indicating that the inequality 

Cj < C2 < 1 for i>3 (24) 

is satisfied for all A^. These results show that, in the critical region of the symmetric phase, 
the two-point Green's function is substantially Gaussian in a large region around fc^, i.e. for 
\k'^ /Mq\ ^ 1, for all A^ from zero to infinity. Another important consequence of the relation 
(P^) is the possibility of evaluating the zero of ^o(y) closest to the origin, ?/o, to a rather 
good approximation by the relationship 

- I/O ^ 1 + C2. (25) 
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The quantity —yo in turn is the scahng hmit of the ratio of the correlation length defined in 
Eq. d^) with the "true" correlation length obtained from the damping factor in the exponen- 
tial long-distance behavior of G{x). Direct calculations of this ratio confirm the relationship 
(^). Similar results have been obtained in the two-dimensional O(A^) models @J^. 

In general models defined on non-Bravais lattices such as the diamond lattice are not 
parity-invariant, and odd-rank operators are allowed in the corresponding expansion of the 
effective Hamiltonian. We finally discuss how space-parity violating terms, when they exist, 
vanish approaching the rotationally- invariant fixed point. This fact should be described by 
a critical exponent pp which should be universal in systems breaking parity at a microscopic 
level, such as ferromagnetic materials having the structure of a diamond lattice. Pp can 
be evaluated on the nearest-neighbor formulation of 0{N) models on the diamond lattice, 
which is not parity invariant. In the correspoding Gaussian theory, or the large- limit of 
O(A^) models, one has Pp = 3, indeed one finds that 

y2xvz G(x,y, z) , . 

'r^, ^ , — . const. 26 

in the critical limit. In general, for finite A^, pp may differ from its Gaussian value. The 
strong-coupling analysis of the odd moments of G{x) on the diamond lattice shows that the 
correction to the Gaussian value of pp is very small. We estimated that < pp — 3 ^ 0.01 
for all > 0. 
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TABLES 



TABLE I. For various values of N, we report the estimates of a obtained by our strong-coupling 

analysis, and by the first nontrivial terms in the 1/N expansion, e-expansion, and (jr-expansion. The 
errors diplayed in the strong-coupling estimates are a rough estimate of the uncertainty, which takes 
into account all the analyses we have performed. 



N s.c. 1/N cxp e-exp 9-exp 

"O 0.00(1) 0.0109 0.0134 

1 0.01(1) 0.0130 0.0157 

2 0.02(1) 0.0140 0.0168 

3 0.03(2) 0.0515 0.0145 0.0170 

4 0.03(2) 0.0386 0.0147 0.0165 
8 0.02(1) 0.0193 0.0137 0.0116 
16 0.009(3) 0.0096 0.0109 0.0096 
32 0.004(2) 0.0048 0.0074 0.0054 
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TABLE II. Strong-coupling estimates of the coefficients C2 and C3 of tlie expansion of go{y), 
and the coefficient di of the expansion of 54 (y). For comparison, we also report estimates from 
the first nontrivial terms in the e-expansion, g'-expansion, and expansion. The errors diplayed 
in the strong-coupling estimates should give an idea of the spread of the results from the various 



Fade— type and integral approximants we considered in our analysis. 
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